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Abstract. I summarize three recent calculations of electromagnetic reactions on deuterium in chiral 
perturbation theory. All of these calculations were carried out to 0(Q 4 ), i.e. next-to-next-to-leading 
order. The reactions discussed here are: elastic electron-deuteron scattering, Compton scattering on 
deuterium, and the photoproduction of neutral pions from deuterium at threshold. 



INTRODUCTION 

Effective field theory (EFT) is a technique commonly used in particle physics to deal 
with problems involving widely- separated energy scales. It facilitates the systematic 
separation of the effects of high-energy physics from those of low-energy physics. In 
strong-interaction physics the low-energy effective theory is chiral perturbation theory 
(%PT) [1]. Here the low-energy physics is that of nucleons and pions interacting with 
each other in a way that respects the spontaneously-broken approximate chiral symmetry 
of QCD. Higher-energy effects of QCD appear in %PT as non-renormalizable contact 
operators. The EFT yields amplitudes which can be thought of as expansions in the ratio 
of nucleon or probe momenta (denoted here by p and q) and the pion mass to the scale 
of chiral-symmetry breaking, A x sb> which is of order the mass of the p meson. The 
existence of a small parameter Q = p/A^sB, g/A^sB, ^jt/A^sb, means that hadronic 
processes can be computed in a controlled way. 

The momentum scale of binding in light nuclei is of order m n and so we should be 
able to calculate the response of such nuclei to low-energy probes using %PT. The result 
is a systematically-improvable, model-independent description. Here I will describe a 
few recent calculations using this approach. Section 2 outlines the %PT expansion, 
and sketches the implications of %PT for calculations of the AW interaction. Section 
3 then looks at electron-deuteron scattering in %PT as a probe of deuteron structre. 
Section 4 examines the use of Compton scattering on the deuteron as a way to extract 
neutron polaraizabilities. I conclude in Section 5 with a discussion of neutral pion 
photoproduction on deuterium at threshold. In these last two reactions I will argue that 
%PT's successful reproduction of the experimental data is in no small part due to its 
consistent treatment of the chiral structure of the nucleon and the deuteron. 



POWER COUNTING AND DEUTERON WAVE FUNCTIONS 



Power counting 

Consider an elastic scattering process on the deuteron whose amplitude we wish to 
compute. If O is the transition operator for this process then the amplitude in question is 
simply with |\|/) the deuteron wave function. In this section, we follow Wein- 

berg [2, 3, 4], and divide the formulation of a systematic expansion for this amplitude 
into two parts: the expansion for O, and the construction of |\|/). 

Chiral perturbation theory gives a systematic expansion for O of the form 

oo 
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where we have labeled the contributions to O by their order n in the small parameter Q 
defined above. Eq. (1) is an operator statement, and the nucleon momentum operator p 
appears on the right-hand side. However, the only quantities which ultimately affect 
observables are expectation values such as For light nuclei this number is 

generically small compared to A x sb- 

To construct one first writes down the vertices appearing in the chiral Lagrangian 
up to order n. One then draws all of the two-body, two-nucleon-irreducible, Feynman 
graphs for the process of interest which are of chiral order Q n . The rules for calculating 
the chiral order of a particular graph are: 

• Each nucleon propagator scales like l/Q; 

• Each loop contributes Q 4 ; 

• Graphs in which both particles participate in the reaction acquire a factor of Q 3 ; 

• Each pion propagator scales like l/Q 2 ; 

• Each vertex from the nth-order piece of the chiral Lagrangian contributes Q n . 

In this way we see that more complicated graphs, involving two-body mechanisms, 
and/or higher-order vertices, and/or more loops, are suppressed by powers of Q. 



Deuteron wave functions 

There remains the problem of constructing a deuteron wave function which is con- 
sistent with the operator O. Weinberg's proposal was to construct a %PT expansion in 
Eq. (1) for the AW potential V, and then solve the Schrodinger equation to find the 
deuteron (or other nuclear) wave function [2, 3, 4]. Recent calculations have shown that 
the NN phase shifts can be understood, and deuteron bound- state static properties reli- 
ably computed, with wave functions derived from %PT in this way [5, 6, 7, 8, 9]. 

Now for %PT in the Goldstone-boson and single-nucleon sector loop effects are 
generically suppressed by powers of the small parameter Q. In zero and one-nucleon 
reactions the power counting in Q applies to the amplitude, and not to the two-particle 
potential. However, the existence of the deuteron tells us immediately that a power 



counting in which loop effects are suppressed cannot be correct for the two-nucleon 
case, since if it were there would be no AW bound state. Weinberg's proposal to instead 
power-count the potential is one response to this dilemma. However, its consistency has 
been vigorously debated in the literature (see [10, 11] for reviews). Recently Beane et 
al. [12] have resolved this discussion, by showing that Weinberg's proposal is consistent 
in the 3 5i — 3 D[ channel. 

One way to understand the %PT power-counting for deuteron wave functions is to 
examine the deuteron wave function in three different regions. Firstly, in the region R 
l/m n the deuteron wave function is described solely by the asymptotic normalizations 
As, A D , and the binding energy B. These quantities are observables, in the sense that they 
can be extracted from phase shifts by an analytic continuation to the deuteron pole. 

The second region corresponds to R ~ 1 /m n . Here pion exchanges play a role in deter- 
mining the AW potential V, and, associatedly, the deuteron wave functions u and w. The 
leading effect comes from iterated one-pion exchange — as has been known for at least 
fifty years. Calculations with one-pion exchange (OPE) defining the potential in this 
regime will be referred to below as "leading-order" (LO) calculations for the deuteron 
wave function. Corrections at these distances come from two-pion exchange, and these 
corrections can be consistently calculated in %PT. They are suppressed by powers of the 
small parameter Q, and in fact the "leading" two-pion exchange is suppressed by Q 2 

relative to OPE. This two-pion exchange can be calculated from vertices in L^J and 
its inclusion in the AW potential results in the so-called "NLO" calculation described 
in detail in Ref. [7]. Corrections to this two-pion-exchange result from replacing one 

(2) 

of the NLO two-pion-exchange vertices by a vertex from . This results in an addi- 
tional suppression factor of Q, or an overall suppression of <2 3 relative to OPE, and an 
"NNLO" chiral potential [5, 6, 7, 8]. More details on this can be found in the contribu- 
tions of Epelbaum and Timmermans to these proceedings. 

Finally, at short distances, R <C l/m n we cutoff the chiral one and two-pion-exchange 
potentials and put in some short-distance potential whose parameters are arranged so as 
to give the correct deuteron asymptotic properties. 



ELASTIC ELECTRON SCATTERING ON DEUTERIUM 

One quantitative test of this picture of deuteron structure is provided by elastic electron- 
deuteron scattering. We thus turn our attention to the deuteron electromagnetic form 
factors Gc, Gq, and Gm- These are matrix elements of the deuteron current i^, with: 

Gc = 3^ y «l|7°|l> + <0|7» + <-l|7°|-l» 5 (2) 

°q = w d ^hz (W )-^ ! 1 )) (3) 
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where we have labeled these (non-relativistic) deuteron states by the projection of the 
deuteron spin along the direction of the momentum transfer q and T| = |q| 2 /(4Mj). 
Gc, Gq, and Gm are related to the experimentally-measured A, B, and T20 in the usual 



way, with T2Q being primarily sensitive to Gq/Gc and B depending only to Gm- Here 
we will compare calculations of the charge and quadrupole form factor with the recent 
extractions of Gq and Gq from data [13]. 

Both of these form factors involve the zeroth-component of the deuteron four-current 
J°. Here we split 7° into two pieces: a one-body part, and a two-body part. The one- 
body part of J° begins at order Q (where we are counting the proton charge e ~ Q) with 
the impulse approximation diagram calculated with the non-relativistic single-nucleon 
charge operator for strutcutreless nucleons. Corrections to the single-nucleon charge 
operator from relativistic effects and nucleon structure are suppressed by two powers 
of Q, and thus arise at 0(Q 3 ), which is the next-to-leading order (NLO) for Gq and 
Gq. At this order one might also expect meson-exchange current (MEC) contributions, 

such as those shown in Fig. 1. However, all MECs constructed with vertices from L^j 
are isovector, and so the first effect does not occur until N 2 LO, or 0(Q 4 ), where an 

NNny vertex from replaces the upper vertex in the middle graph of Fig. 1, and 
produces an isoscalar contribution to the deuteron charge operator. (This exchange- 
charge contribution was first derived by Riska [14].) 



FIGURE 1. The impulse-approximation contribution to Gc and Gq is shown on the left, while two 
meson-exchange current mechanisms which would contribute were the deuteron not an isoscalar target 
are depicted in the middle and on the right. 

The most important correction that arises at NLO is the inclusion of nucleon structure 
in %PT. At 0(Q i ) the isoscalar form factors are dominated by short-distance physics, and 
so the only correction to the point-like leading-order result comes from the inclusion of 
the nucleon's electric radius, i.e. 

^' ) X PTNLO = 1 -^4 5)2 )^ 2 - (5) 

This description of nucleon structure breaks down at momentum transfers q of order 300 
MeV. There is a concomitant failure in the description of eD scattering data [15, 16]. 
Consequently, in order to focus on deuteron structure, in the results presented below 
I have chosen to circumvent this issue by using a "factorized" inclusion of nucleon 
structure [16]. This facilitates the inclusion of experimentally-measured single-nucleon 
form factors in the calculation, thereby allowing us to test how far the theory is able to 
describe the two-body dynamics that takes place in eD scattering. 

The results for Gc and Gq are shown in Fig. 2. The figure demonstrates that conver- 
gence is quite good below q ~ 700 MeV — especially for Gc- The results shown are for 
the NLO chiral wave function, but the use of the NNLO chiral wave function, or indeed 
of simple wave functions which include only one-pion exchange, do not modify the pic- 
ture greatly below q = 700 MeV [16]. It is also clear that — provided information from 
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FIGURE 2. The deuteron charge and quadrupole form factors to order Q 4 in chiral perturbation theory. 
The experimental data is taken from the compilation of Ref. [13]. Gq is in units of fm 2 . 



eN scattering is taken into account — %PT is perfectly capable of describing the charge 
and quadrupole form factors of deuterium at least as far as the minimum in Gc- This 
result is extremely encouraging for the application of EFT to light nuclei. 

Gm can be obtained in a similar way, but, importantly, the LO contribution to Gm 
is 0(Q 2 ). Furthermore, no two-body mechanism enters until 0(Q 5 ), when an undeter- 
mined two-body counterterm appears [15, 17]. Results for F M at 0(Q 4 ) turn out to be 
of similar quality to those for Fq [15, 16], but are somewhat more sensitive to short- 
distance physics, as was expected given the presence of the 0(Q 5 ) counterterm in this 
observable. 

The static properties of the deuteron obtained in this expansion are also generically 
quite reasonable, and have good convergence properties. The one exception to this is 
the deuteron quadrupole moment, Q c i, which, is underpredicted by about 4% — as is also 
true in all modern potential-model calculations [18]. However, such an underprediction 
is not unexpected in %PT since simple estimates of the effect on Qj of higher-order terms 
in the chiral Lagrangian suggest that a discrepancy of order 5% is to be expected at 
0(Q 4 ). Qd is rather sensitive to short-distance physics and it transpires that higher-order 
counterterms have a larger effect on it than on Gc [19, 20]. This way of understanding 
the "<2j puzzle" is one example of the way in which %PT can assist in the analysis 
of electromagnetic currents for few-body systems. For an analogous application of 
%PT/EFT ideas to the important solar reaction pp — > de + V e see Ref. [21] and the 
contribution of Marcucci to these proceedings. 



COMPTON SCATTERING ON DEUTERIUM 



Compton scattering on the nucleon at low energies is a fundamental probe of the 
long-distance structure of these hadrons. This process has been studied in %PT in 



Ref. [22, 23], where the following results for the proton polarizabilities were obtained 
atLO: 



h = Jffh = 12 - 2 x KT 4 fm 3 ; p\ = n jY* = 1.2 x lO^fm 3 . 
p 384tu 2 ./>* Fp 76Sn 2 f%m n 



Recent experimental values for the proton polarizabilities are [24] 
a p + $ p = 13.23 ± 0.86l[j;2g x 10~ 4 fm 3 , 

a p - |3 p = 10.11 ± 1.741S x 10- 4 fm 3 , (7) 

where the first error is a combined statistical and systematic error, and the second set of 
errors comes from the theoretical model employed. These values are in good agreement 
with the %PT predictions. 

Chiral perturbation theory also predicts a n = a p , p\ 3 = (3 p at this order. The neutron 
polarizabilities a n and (3„ are difficult to measure, due to the absence of suitable neutron 
targets, and so this prediction is not well tested. One way to extract a n and |3 n is 
to perform Compton scattering on nuclear targets. Coherent Compton scattering on a 
deuteron target has been measured at Ey = 49 and 69 MeV by the Illinois group [25] and 
Ey = 84.2 — 104.5 MeV at Saskatoon [26]. The amplitude for Compton scattering on the 
deuteron clearly involves mechanisms other than Compton scattering on the individual 
constituent nucleons. Hence, the desire to extract neutron polarizabilities argues for a 
theoretical calculation of Compton scattering on the deuteron that is under control in the 
sense that it accounts for all mechanisms to a given order in %PT. 




FIGURE 3. Graphs which contribute to Compton scattering on the deuteron at o(Q 2 ) (a) and o(Q i ) 
(b-d). The sliced and diced blobs are from L^j (c) and (d). Crossed graphs are not shown. 

The Compton amplitude we wish to evaluate is (in the yd center-of-mass frame): 

where M (M 1 ) is the initial (final) deuteron spin state, and X (X!) is the initial (final) 
photon polarization state, and k (k 1 ) the initial (final) photon three-momentum, which 
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are constrained to \k\ = \k \ = (£>. The amplitude T yN represents the graphs of Fig. 3 and 
Fig. 4b where the photon interacts with only one nucleon. The amplitude T 2 ^ N represents 
the graphs of Fig. 4a where there is an exchanged pion between the two nucleons. 
The LO contribution to Compton scattering on the deuteron is shown in Fig. 3(a). 

This graph involves a vertex from and so is 0(Q 2 ). This contribution is simply the 
Thomson term for scattering on the proton. There is thus no sensitivity to either two- 
body contributions or nucleon polarizabilities at this order. At 0(Q i ) there are several 
more graphs with a spectator nucleon (Figs. 3(b),(c),(d)), as well as graphs involving 
an exchanged pion with leading order vertices (Fig. 4(a)) and one loop graphs with a 
spectator nucleon (Fig. 4(b)) [27]. Graphs such as Fig. 4(b) contain the physics of the 
proton and neutron polarizabilities at 0(Q 3 ) in %PT. 
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FIGURE 4. Graphs which contribute to Compton scattering on the deuteron at (Q 3 )- Crossed graphs 
are not shown. 

We employed a variety of wave functions \|/, and found only moderate wave-function 
sensitivity. Results shown here are generated with the NLO chiral wave function of 
Ref. [7]. Fig. 5 shows the results at Ey — 49, 69, and 95 MeV. For comparison we have 
included the calculation at 0(Q 2 ) in the kernel, where the second contribution in Eq. (8) 
is zero, and the single-scattering contribution is given solely by Fig. 3(a). At 0(Q 3 ) all 
contributions to the kernel are fixed in terms of known pion and nucleon parameters, 
so to this order %PT makes predictions for deuteron Compton scattering. We also show 
the 0(Q 4 ) result which will be discussed below. The curves indicate that higher-order 
corrections get larger as CO is increased — as expected. 
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FIGURE 5. Results of the 0(Q 2 ) (dotted line), 0{Q 3 ) (dashed line), and 0(Q 4 ) (solid line) calculations 
for E y = 49 MeV, 69 MeV and 95 MeV respectively from left to right. 

We have also shown the six Illinois data points at 49 and 69 MeV [25] and the 
Saskatoon data at 95 MeV [26]. Statistical and systematic errors have been added in 
quadrature. It is quite remarkable how well the 0(Q 2 ) calculation reproduces the 49 



MeV data. However, the agreement is somewhat fortuitous: there are significant 0(Q 3 ) 
corrections. Note that at these lower photon energies Weinberg power counting begins to 
break down, since it is designed for (0 ~ %, and does no? recover the deuteron Thomson 
amplitude as (0 — > 0. Correcting the power counting to remedy this difficulty appears to 
improve the description of the 49 MeV data, without significantly modifying the higher- 
energy results [27]. Meanwhile, the agreement of the 0(Q 3 ) calculation with the 69 
MeV data is very good, although only limited conclusions can be drawn. These results 
are not very different from other, potential-model, calculations [28, 29, 30]. They are 
also quite similar to those obtained in AW EFTs without explicit pions (see Ref. [31], 
and the contribution of GrieBhammer to these proceedings). However our calculation is 
the only one that does not employ the polarizability approximation for the yN amplitude. 

At 0(Q 4 ) single-nucleon counterterms which shift the polarizabilities enter the cal- 
culation. However, there are still no two-body counterterms contributing to yd — > yd at 
this order. In this sense Compton scattering on deuterium at 0(Q 4 ) is analogous to the 
reaction yd — > n°d discussed below: an 0(Q 4 ) calculation allows us to test the single- 
nucleon physics which is used to predict the results of coherent scattering on deuterium, 
since there are no undetermined parameters in the two-body mechanisms that enter to 
this order in the chiral expansion. 

The 0(Q 4 ) plot shown above is a partial calculation at that order. It includes all two- 
body mechanisms at 0(Q 4 ), and some one-body mechanisms, the latter being calculated 
in ways motivated by dispersion-relation analyses. The values of the proton polarizabil- 
ities used in Fig. 5 were taken from Eq. (7). Meanwhile the (disputed) neutron-atom 
scattering value of Ref. [32] was employed for a n and the Baldin sum rule used to fix 
pV To demonstrate the sensitivity to a n , we can also use the freedom in the single- 
nucleon amplitude at 0(Q 4 ) to fit the SAL data using the incomplete 0(Q 4 ) calculation. 
A reasonable fit at backward angles can be achieved with a n = 4.4 and % = 10. These 
numbers are in startling disagreement with the 0(Q 3 ) %PT expectations. We have also 
plotted the cross-section at 69 MeV with a n = 4.4 and |3„ = 10: this curve misses the 
Illinois data. (Similar results were found in a potential model in Ref. [29].) This sit- 
uation poses an interesting theoretical puzzle. A full 0(Q 4 ) calculation in yPT using 
the recently derived single-nucleon Compton amplitude [33] is necessary before firm 
conclusions can be drawn. Such a calculation is in progress [34]. 
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FIGURE 6. Results of 0(Q 3 ) (dashed line), partial 0(Q 4 ) (solid line), and partial 0(Q 4 ) with modified 
a„ (dot-dashed line), xPT calculations for E y = 69 MeV (left panel) and 95 MeV (right panel). 



NEUTRAL PION PHOTOPRODUCTION ON DEUTERIUM 



Pion photoproduction on the nucleon near threshold has been studied up to 0(Q 4 ) 
in %PT with the delta integrated out in Ref. [35]. The differential cross-section at 
threshold is given solely by Eq + , the electric dipole amplitude. In %PT neutral pion 
photoproduction on the nucleon does not begin until 0(Q 2 ), where there is a tree- 
level contribution from a jkNN vertex. Then at 0(Q 3 ) there are tree-level contributions 
involving the magnetic moment of the nucleon. The sum of these 0(Q 2 ) and 0(Q 3 ) 
effects reproduces an old "low-energy theorem". Also at 0(Q 3 ) there occur finite loop 
corrections where the photon interacts with a virtual pion which then rescatters on the 
nucleon. This large quantum effect was missed in the old "low-energy theorem" and is 
absent in most models. At 0(Q 4 ) there are loops with relativistic corrections, together 
with a counterterm. The proton amplitude that results from fitting this counterterm to 
data produces energy-dependence in relatively good agreement with the recent results 
from Mainz and Saskatoon. At the same order there is also a prediction for the near- 
threshold behavior of the yn — > n°n reaction; the cross-section is considerably larger than 
that obtained in models that omit the important pion-cloud 0(Q 3 ) diagrams. In fact, the 
0(Q 4 ) prediction is that the cross section for yn — > n°n is about four times larger than 
that for neutral pion photoproduction on the proton. A good way to test this prediction 
is to study neutral pion photoproduction on the deuteron. 

If deuterium is to be used as a target in this particular process then two-body mech- 
anisms that contribute to the reaction must be calculated. %PT provides an ideal way to 
do this, as was demonstrated by Beane and collaborators [36]. 

Consider the reaction y(k) +d(p\) — > Tt°(q) + d(p2) in the threshold region, q ~ 0, 
where the pion is in an S wave with respect to the center-of-mass (cm) frame. For real 
photons the threshold differential cross section is: 



\k\ do 
\q\ dQ. 



l?l=o 



= \E 2 d . (9) 



We now present the results of the chiral expansion of the dipole amplitude E d to 0(Q 4 ). 
Two-body contributions to E d do not begin until 0(Q 3 ). Thus, to 0(Q 4 ) we have 

E d = (¥l^l¥) + (¥l^ 3 V) + +(¥l^V) (10) 
= E s d s + E d b > 3 +E d b > 4 , (11) 

where we have explicitly isolated the 0(Q 3 ) and 0(Q 4 ) two-body contributions to the 
operator O, and |\|/) is a deuteron wave function. 

The single- scattering contribution to E d , E d , is given by all diagrams where the 
photon is absorbed and the pion emitted from one nucleon with the second nucleon 
acting as a spectator, i. e. the impulse approximation result. Since the %PT results for 
the elementary 5-wave pion production amplitudes to 0(Q 4 ) [35] are known there is an 
0(Q 4 ) %PT prediction for E d . This is the prediction we want to test. 

But, before comparing this prediction with the experimental data we must compute 
the two-body contributions to E d . Those of 0(Q 3 ) which survive at threshold, in the 
Coulomb gauge, are shown in Fig. 7 [37]. 




(a) 
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FIGURE 7. Two-nucleon graphs which contribute to neutral pion photoproduction at threshold to 
0(Q i ) (in the Coulomb gauge). All vertices come from L^j. 



At 0(Q 4 ), we have to consider the two-nucleon diagrams — some of which are shown 

Fig. 8, with the blob characteriz 
corrections to the graphs in Fig. 7. 



(2) 

in Fig. 8, with the blob characterizing an insertion from L^. There are also relativistic 




FIGURE 8. Characteristic two-nucleon graphs contributing at 0(Q A ) to neutral pion photoproduction. 



The hatched circles denote an insertion from . 



One can show that the only terms that survive at threshold result in insertions ~ 1/2M, 
~ ^a/2M and ~ Ko.i. To 0(Q 4 ) there are no four-nucleon operators contributing to the 
deuteron electric dipole amplitude, and so no new, undetermined parameters appear. 
The only free parameter is fixed in yp — > K p, and so a genuine prediction can be, and 
was, made for the reaction yd — > n°d at threshold. 

We now present results. Using a variety of deuteron wave functions the single- 
scattering contribution is found to be [36]: 

E s d s = (0.36 ±0.05) x \Q- 3 /m n+ . (12) 

The sensitivity of the single-scattering contribution E S J to the elementary neutron am- 
plitude may be parameterized by: 



T?SS 



0.36-0.38-(2.13-£^' J ) x lO^/m^. (13) 



The two-nucleon contribution is evaluated at 0(Q 3 ) using different deuteron wave 
functions [36]. The results prove to be largely insensitive to the choice of wave function, 
as do the two-body contributions at 0(Q 4 ). Choosing the AV18 wave function for 



definiteness the results are summarized in Table l 1 . The 0(Q 4 ) contributions give 
corrections of order 15% to the 0(Q 3 ) two-nucleon terms. We also observe that E'J'' 4 
is of the same size as E s /, clearly demonstrating the need to go to this order in the 
expansion. This gives us confidence that the %PT expansion is controlled, and that we 
may compare experimental data with the 0(Q 4 ) prediction [36] : 

Ed = (-1.8 ±0.2) x 10~ 3 /m^. (14) 



TABLE 1. Values for Ed in units of 10 3 /m Jt + from one- 
nucleon contributions (IN) up to 0(Q 4 ), two-nucleon kernel 
(2A0 at 0(Q 3 ) and at 0(Q 3 ) , and their sum (IN + 2N). 



IN 


2N 


1N + 2N 


Q + Q 2 + Q 3 + Q 4 


q' e 4 


Q + Q 2 + Q 3 + Q 4 


0.36 


-1.90 -0.25 


-1.79 



To see the sensitivity to the elementary neutron amplitude, we set the latter to zero and 
find E d = —2.6 x 10~ 3 /m % + (for the AV18 potential). Thus %PT makes a prediction that 
differs markedly from conventional models, with that difference arising predominantly 
from a different result for E^. An experimental test of this prediction was carried 
out recently at Saskatoon [38]. The results for the pion photoproduction cross-section 
near threshold are shown in Fig. 9, together with the %PT prediction at threshold. The 
agreement with %PT to order 0(Q 4 ) is not better than a reasonable estimate of higher- 
order terms, but %PT is clearly superior to models. This is compelling evidence of the 
importance of chiral loops, and also testimony to the consistency and usefulness of %PT 
in analyzing low-energy reactions on deuterium. 
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FIGURE 9. Reduced cross-section Or = (k/q)o in /jb for neutral pion photoproduction as function of 
the photon energy in MeV. Threshold is marked by a dotted line. Squares are data points from Ref. [38] 
and the star is the %PT prediction of Ref. [36]. Figure courtesy of Ulf MeiBner. 



Similar answers were obtained with the wave functions of Ref. [5], which should be equivalent to the 
NNLO wave function discussed above, up to higher-order terms. 
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